Rapid Suppression of the Spin Gap in Zn-doped CuGeC>3 and SrCu 2 03 
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The influence of non-magnetic impurities on the spectrum and dynamical spin structure factor of 
a model for CuGeC>3 is studied. A simple extension to Zn-doped SrCu2 03 is also discussed. Using 
Exact Diagonalization techniques and intuitive arguments we show that Zn-doping introduces states 
in the Spin-Peierls gap of CuGeC>3. This effect can be understood easily in the large dimerization 
limit where doping by Zn creates "loose" S=l/2 spins, which interact with each other through 
very weak effective antiferromagnetic couplings. When the dimerization is small, a similar effect 
is observed but now with the free S=l/2 spins being the resulting S=l/2 ground state of severed 
chains with an odd number of sites. Experimental consequences of these results are discussed. It 
is interesting to observe that the spin correlations along the chains are enhanced by Zn-doping 
according to the numerical data presented here. As recent numerical calculations have shown, 
similar arguments apply to ladders with non-magnetic impurities simply replacing the tendency to 
dimerization in CuGeOa by the tendency to form spin-singlets along the rungs in SrCu2 03. 

PACS numbers: 64.70.Kb,75.10.Jm,75.50.Ee 



I. INTRODUCTION 

The recent discovery of a spin-Peierls (SP) transition 
in the inorganic compound CuGe03 has received consid- 
erable attention. The effect is presumably caused by 
the coupling of the spins along the Cu — O2 chains with 
three-dimensional phonons. This coupling dimerizes the 
chains reducing (increasing) the lattice spacing for, e.g., 
the even (odd) links, and thus promoting the formation of 
spin singlets along the short bonds. A spin-gap is formed 
in the spectrum, which is roughly the energy necessary 
to break the even-bond spin singlets. Soon after the dis- 
covery of this material, studies where Cu 2+ (S=l/2) was 
substituted by Zn 2+ , Ni 2+ or Mn 2+ were reported. Q 
Based on magnetic susceptibility % data analysis, a rapid 
suppression of the SP spin-gap was observed, with a small 
2% concentration of Zn being sufficient to induce a transi- 
tion to a new phase, x, specific heat and inelastic neutron 
scattering (INS) measurements ||-|| have unambiguously 
shown that this new phase is a three dimensional (3D) 
Neel state. The simultaneous existence in the experimen- 
tal data of features indicating SP and Neel order has also 
been noticed. |3| Similar results have been reported with 
Si substituting Ge, § and Ni,Mn,Mg replacing Cu. §@] 
These results are apparently inconsistent with the gapless 
SP state predicted by mean- field theory. § Recently @, 
INS measurements in good quality samples of the tem- 
perature dependence of the superlattice peaks caused by 
the lattice dimerization were interpreted as evidence that 
the SP actually survives up to 6% Zn-doping. 

The experimental evidence accumulated recently in 
doped CuGe03 shows two interesting phenomena namely 



(i) the unexpectedly rapid suppression of the SP order 
by Zn-doping, and (ii) its replacement by 3D AF or- 
der. Since INS results for 2% Zn-doping ||] show already 
an important reduction of the SP transition at temper- 
atures well above those where AF order starts, we con- 
jecture that these two phenomena could be studied in- 
dependently. It is remarkable that similar experimental 
results jl(]] have been observed in a quasi one-dimensional 
"ladder" system, Sr(Cui_ x Zn x ) 2 03, where a spin-gap ex- 
ists at x=0 due to the formation of spin-singlets along 
the ladder rungs, [|ll| without lattice deformations. This 
suggests that the rapid suppression of a spin-gap by Zn- 
doping and its subsequent replacement by 3D Neel order 
may be a general phenomenon independent of the origin 
of such a spin-gap. 

In this paper, the effect of non-magnetic Zn-doping on 
the SP state of CuGe03 is studied. A simple general- 
ization to the case of Zn-doped ladders is also discussed. 
Using numerical techniques we show that the SP spin- 
gap is rapidly suppressed with Zn-doping, as observed 
experimentally. The effect is caused by unpaired spins 
created by the breaking of spin-singlets upon doping (in 
agreement with conclusions reached in recent numerical 
studies of even- leg ladder systems p2|). These unpaired 
spins coupled forming an effective random spin- 1/2 chain. 
Our analysis does not include intra-chain couplings that 
could stabilize the 3D AF order, analysis that is post- 
poned for a future publication [[13[ (however, in Sec. VII 
below we discuss the enhancement observed numerically 
of spin correlations upon Zn doping. This effect may 
contribute to the stabilization of the 3D AF order in real 
materials). To study the dimerization of a spin-1/2 chain 
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it is natural to use the standard Heisenberg model includ- 
ing a static modulation of the nearest-neighbor (NN) ex- 
change to account for the lattice distortion. In addition, 
recent studies have shown that a next-nearest-neighbor 
(NNN) spin-spin interaction is also needed to properly 
fit experimental susceptibility data. |l^-|l7|] This is rea- 
sonable based on the structure of CuGeC>3 where Cu-Cu 
interactions bridged by two oxygens, Cu-O-O-Cu, may 
not be negligible. Thus, the model used here is 

TL = J\ ^ Sj • Sj + J 2 ^ Sj • Sfe — 

(ij) ({ik)) 

J 1( ^(-irS 4 -S,, (1) 

(ij) 

where (ij) (((ik)}) denote NN (NNN) sites, S regulates 
the amount of dimerization, and J\ , Ji are the NN and 
NNN Heisenberg couplings, respectively. The rest of the 
notation is standard. A modulation of Ji is a higher or- 
der effect neglected here. To model CuGeC>3, we here 
use 6 = 0.03 and J2/J1 = 0.24, |f| but other set of pa- 
rameters could be chosen. |l4|] The qualitative results de- 
scribed below do not depend crucially on the actual val- 
ues of 6 and J2/ J\. The introduction of a non-magnetic 
impurity at site i is mimicked in Eq.(l) by removing the 
Ji, Ji and (5-terms associated with this site. Note that 
the Ji NNN coupling between sites i—1 and i + 1 is not 
influenced by an eventual impurity located at i. This cou- 
pling across impurities plays an important role in linking 
chain segments severed by Zn-doping. 



II. LARGE DIMERIZATION LIMIT 

To guide the intuition on how Zn-doping affects the 
properties of model Eq.(l), it is convenient to study the 
limit where the dimerization 6 is large, inducing a well- 
defined pattern of dimers (spin-singlets) in the ground 
state of the system. To illustrate the physics of the prob- 
lem consider in Fig.l a simple case where two impurities 
are randomly distributed on an 8-site chain using peri- 
odic boundary conditions (PBC). Impurity- 1 is located 
at site 1, and the second one at different sites along the 
chain. There are six independent cases shown in Fig.l. 
In Fig. 2a, the longitudinal dynamical spin structure fac- 
tor S zz (tt, ui), which is the spectral decomposition of the 
operator (1/N) Y^i(~^-T^i > ^ s shown for the different po- 
sitions of the impurities (N is the number of sites). We 
consider momentum 7r because in this subspace is where 
the smallest gap is found in the undoped system. The 
calculation was performed with standard Exact Diago- 
nalization (ED) techniques, Q and using <5 = 0.5 (pro- 
viding a strong exchange J = 1.5 Jj for half the links of 
the chain, and a weak exchange J — 0.5 Ji for the other 
half). We observed that for such a large dimerization 
finite size effects are small. We also verified that using 



chains with open boundary conditions (OBC) the con- 
clusions of this discussion are unchanged. 
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FIG. 1. Eight site chain with PBC and two impurities 
located in the positions of the full dots. The thick (thin) 
links denote a large (small) spin-spin coupling of magnitude 
Ji(l + S) (Ji(l — S)). The dashed line linking some pairs of 
sites illustrates the presence of the J2 coupling across impu- 
rities. The other J2 couplings are not shown explicitly. The 
case of impurities at sites 1 and 7 is identical to having im- 
purities at 1 and 3. 

Let us consider the six impurity distributions depicted 
in Fig.l: (1) for impurities at sites 1-2 a strong dimer 
is entirely removed. All the other spins are themselves 
strongly coupled in dimers, and thus we expect spin exci- 
tations at energies u> ~ 1.5Ji which are needed to trans- 
form a spin singlet into a triplet, in good agreement with 
Fig. 2a. In all the other cases below, weight at this large 
uj is also expected due to the presence of strong dimers in 
the ground state; (2) if the impurities are at sites 1 and 
3, they break two dimers producing two unpaired spins. 
These spins tend to form a spin-singlet using the coupling 
Ji , since the rest of the spins are mostly forming dimers 
as shown in Fig.l. Then, singlet-triplet excitations are 
expected at u ~ 0.24Ji. The anomalous small weight 
of the associated pole in Fig. 2a is caused by the use of 
momentum 7r in the dynamical spin structure factor that 
gives small weight to singlets where both spins are at even 
or odd sites; (3) with impurities at 1 and 4, two strong 
dimers are erased, and the two remaining spins are linked 
into a singlet through J = 0.5Ji inducing a pole at en- 
ergy co ~ 0.5 Ji as in Fig. 2a; (4) when the impurities are 
at opposite sites of the chain, two free spins are created 
which are linked by a very weak effective exchange since 
two of the intermediate spins form a strong singlet. This 
would imply the presence of weight at small u>. However, 
sites 2 and 6 are both even and the spin structure fac- 
tor at q = 7r suppresses the signal; (5) with impurities at 
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sites 1-6, the two free spins are also weakly coupled into a 
singlet (weight at small w is thus expected in agreement 
with Fig. 2a), but now being at even and odd sites their 
contribution to S zz (tt,uj) is not suppressed by the sign 
modulation. The effective coupling is J e // ~ 0.07 Ji; (6) 
finally, for impurities at 1 and 8, the free spins are again 
weakly linked and the intensity of the pole signaling their 
singlet-triplet excitation is enhanced by the momentum 
q = n. Thus, weight is expected at small lo, as in Fig. 2a. 
Note that in this case the effective exchange is the small- 
est of the cases shown in Fig.l, J e ff ~ 0.02 J 1; since it is 
mediated by two strong dimers. 
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FIG. 2. (a) S zz (n,ui) vs ui on an eight site chain with PBC 
and using a dimerization 8 — 0.5. The impurities are dis- 
tributed as in Fig.l. The 5-functions have been given a width 
e — 0.03. The labels indicate the position of the second im- 
purity, with the first one located at site 1; (b) S^S (t, w) (av- 
erage of results Fig. 2a over the positions of the impurities) 
vs lu (solid line). The dashed line indicates the result in the 
absence of impurities. 

In Fig. 2b, the equal- weight average of S zz (ir,ui) over 
the impurity positions, denoted as S zz (it,uj), is shown. 
The result in the absence of impurities is also presented. 
It is clear that the introduction of impurities broadens the 
main peak, and, more importantly, it introduces weight 
inside the original spin-gap. The peak the smallest in 
intensity is caused by the special case when two impuri- 
ties are at close distance, while the weight at very small 
lo is originated by the weak coupling into singlets of un- 
paired spins at large distances. The weight inside the gap 
should diminish as the density of impurities is reduced. 
Note that the effective interaction between spins at large 
distances can be understood using a decimation proce- 
dure analog to that employed by Dasgupta and Ma [jl9| , 
in their study of random spin-1/2 chains. 



III. REALISTIC PARAMETERS 

In Fig. 3a, S zz (tt,lo) is shown as 5 varies from 0.5 to 
0.03 now using a chain of 20 sites and still 2 impurities. 



The peak at high to moves down in energy as 6 is reduced, 
as expected since it is originated by the singlet-triplet ex- 
citations of the dimers. On the other hand, the weight at 
small u> increases slightly its energy as S is reduced. This 
is also compatible with the expectation that this feature 
is created by the AF coupling between "free" spins lo- 
cated at large distances from each other. A smaller S 
increases the coupling between those spins since the ten- 
dency to form tight intermediate spin singlets is reduced. 
The evolution with 6 of Fig. 3a is smooth and suggests 
that the intuitive explanation for the filling of the gap 
found at large 5 could be qualitatively correct even in 
the more realistic regime of 6 <C 1. 
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FIG. 3. (a) S%%(n,(j) vs lo on a 20-site chain parametric 
with 8 as indicated. The peaks in the spectral weight have 
been given a width e = 0.03. (b) £^(7^0;) vs u for 8 — 0.03 
with a resolution e = 0.006. The labels of the peaks refer 
to the position of a second impurity on a chain of 20-sites 
with the first impurity located at site 1. The missing sites are 
related to those shown by symmetry, or they have negligible 
weight. 

On Fig. 3b, S zz (tt,lo) is shown for 5 — 0.03 but with a 
higher resolution than used in Fig. 3a. Each peak repre- 
sents the position of an excited state of the system with 
spin 1. The labels refer to the position of the second 
impurity, the first one being always at site 1. With this 
resolution, we can identify the origin of the two dominant 
features observed in the spin structure factor Fig. 3a at 
small S, namely the presence of low-w and high-w spec- 
tral weight. Fig. 3b shows a clear difference between even 
and odd sites. To understand this effect, note that the 
case where the second impurity is on an odd-site corre- 
sponds to a partition of the chain into two segments each 
with an odd number of spins. This effectively produces 
two S=l/2 states in each segment which are coupled by 
the relatively small J2. The fact that such an effective 
"free" spin-1/2 in each segment is now spread over more 
than one lattice spacing contributes to the reduction of 
the effective exchange leading to the very small energy 
excitation of, e.g., the case where the impurities are at 



3 



sites 1 and 9. Thus, the same reasoning leading to the 
appearance of weight at small u> in the large 5 limit, can 
be repeated for small S still using J 2 as the link between 
S=l/2 effective spins. 



IV. INTUITIVE PICTURE 

The basic idea found numerically in the previous sec- 
tion is that Zn-doping divides a S=l/2 chain into seg- 
ments each carrying either a total spin 1/2 or 0, depend- 
ing on the number of spins in the segment. Then, at low 
temperatures each segment can be replaced by just one 
effective spin with coupling that depend on the distance 
between impurities (i.e. on the length of the segments). 
This is reminiscent of a similar approach followed recently 
to study spin-1/2 chains with randomly distributed fer- 
romagnetic and antifcrromagnetic bonds, as realized in 
Sr3CuPto.5Iro.5O6- p3fl If such a simple picture were cor- 
rect for Zn-doped CuGe03, then the low-energy states of 
the complete doped chain could be calculated using an ef- 
fective Heisenberg Hamiltonian that couples the effective 
neighboring spins. 
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FIG. 4. (a,b,c) Clusters used to illustrate the main intuitive 
ideas explained in the text. Each of the three segments carry 
an effective spin 1/2, denoted by Si, S2 and S3. The circles 
are the impurities. The thick (thin) bonds have large (small) 
exchanges. The dashed line denotes some of the J2 couplings. 

Let us verify this scenario with the example of a 15- 
sites chain with OBC, and two impurities located at sites 
4 and 10 dividing the lattice into three segments, each 
with an odd number of spins (Fig. 4a). This setup is sup- 
posed to produce three effective spin-1/2 that we denote 
as Si, 5*2 and S3. Its low energy Hamiltonian should be 

7~i = J12S1 ■ S2 + J23S2 • S3, (2) 
with J12 and J23 denoting effective exchange couplings, 



which in principle correspond to the singlet-triplet energy 
separation Agj> for the case where two adjacent segments 
of Fig. 4a are analyzed in isolation. Diagonalizing exactly 
these coupled segments of 3 and 5 sites (Fig. 4b) and 5 
and 5 sites (Fig. 4c), respectively, J12 and J23 have been 
calculated, and the resulting couplings used in the full 
three-segment 15-site problem of Fig. 4a, which can also 
be exactly diagonalized to gauge the accuracy of Eq.(2). 
Fig. 5 shows that the agreement between the exact result 
and the effective model is excellent providing support to 
the simple picture outlined above. 
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FIG. 5. Energy levels of the effective Hamiltonian Eq.(2) 
compared against those of the 15-site cluster shown in Fig. 4a. 
The eigenvalues of Eq.(2) were shifted by a constant such 
that its ground state energy coincides with that of the 15-site 
chain. 

To complete this analysis, we must investigate the case 
where the segments have an even number of spins. Here 
in the large 6 limit there are two clearly different sit- 
uations: 1. all the spins of the even-site segment are 
coupled into strong dimers or 2. there are two "loose" 
spins at the ends of the segment. In the first case, it 
is clear that an even-site mostly dimerized segment in- 
tercalated between two odd-spin chains simply reduces 
substantially the coupling between the effective spins of 
such odd-chains (with the effect becoming more drastic 
as the length of the even-chain grows). This implies that 
the random distribution of antiferromagnetic couplings 
resulting from this effect may contain large weight at 
small exchanges i.e. it will not be centered at J\ or J2. 
In the second case where there are two loose spins in the 
same segment, their mutual interaction is also weak since 
the rest of the spins form mostly singlets at least at large 
8. However, the free spins at the ends of the segment 
can couple through J2 with the free spins of the adjacent 
chains. We have carried out explicit ED calculations that 
indeed confirm this simple picture. 
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V. COMPARISON WITH AVAILABLE 
EXPERIMENTS 



tested by INS can survive the introduction of doping up 
to large impurity concentrations. However, \ denned as 



In Fig. 6a, the dynamical spin structure factor 
S*v(q, uj), corresponding to the spectral decomposition of 
the operator (1/N) ^ z e lql Sf , is shown parametric with 
momentum q. As discussed before, this is an "equal- 
weight" average over the positions of the impurities. Note 
that in real materials we may have some other probability 
distribution of hole positions determined by energy and 
entropy factors. Keeping this detail in mind, we proceed 
with the analysis of Fig. 6a assuming an equal-weight dis- 
tribution. For all momenta there is a clear distinction 
between the low and high-w features discussed before. 
Fig. 6b contains the g-dependence of the centroid of the 
low and high-w spectral weight. It is clear that the fea- 
tures at large u> are remnants of the S=l excitations of the 
Ji -Ji model (here S = 0.03 plays only a secondary role). 
||16| These results also clearly show that weight is gen- 
erated for energies inside the original SP spin-gap of the 
undoped system upon non-magnetic doping of the chains. 
This weight has a mild q-dependence. Thus, based on the 
intuitive picture developed in this paper the spin-gap is, 
rigorously speaking, destroyed for an infinitesimal con- 
centration of Zn-doping. However, the weight inside the 
gap is also proportional to the impurity concentration 
and it may not be experimentally detectable until a fi- 
nite doping threshold is reached. 
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FIG. 6. (a) £>2J(<2,a>) vs u) for a 20-site chain with 2 im- 
purities, parametric with q. The triangle denotes the rough 
position of the low and high-tj spectral weight; (b) Position of 
the dominant high and low-w features of (a) compared with 
results in the absence of impurities. 

It is important to remark that even though the gap 
is closed with doping, in Fig. 6a clear remnants of the 
original undoped system are observed. Then, accord- 
ing to our picture INS experiments should indicate the 
presence of dimerization and AF correlations for Zn- 
concentrations where other probes, specially magnetic 
susceptibility measurements, indicate the disappearance 
of the spin-gap. Short-distance dimerization correlations 



is dominated by small energy excitations even if they 
carry small weight. Using q = 7r/10 in Eq.(3), we calcu- 
lated numerically \ a t zero temperature using a 20-site 
chain with and 2 impurities. The ratio Xdoped I Xundoped 
is about 40. Therefore this explains why x measurements 
can show no trace of a SP phase, while INS still suggest 
its presence. || It is the short- against long-character of 
the correlations that establish the difference, x reacts to 
the long distance structure, while INS captures fluctua- 
tions at all distances. 



VI. EXPERIMENTAL PREDICTIONS 

The numerical and analytical results discussed in this 
paper suggest that CuGe0 3 looses its gap immediately 
after an infinitesimal amount of Zn replaces Cu. This 
dramatic effect produces interesting experimental conse- 
quences in addition to those described in the previous sec- 
tion. The localized spin- 1/2 states in the vicinity of the 
nonmagnetic ions interact weakly through a Heisenberg 
Hamiltonian with antiferromagnetic couplings having a 
strength that depends on the distance between the Zn 
impurities. Since this distance is random, the couplings 
are random and locate the low energy behavior of Zn- 
doped CuGeC>3 in the universality class of the random 
exchange Heisenberg model. This low energy effective 
model, defined by the Hamiltonian 

H = 2^ Jij Si ' Sj 

in a standard notation, has been analyzed using decima- 
tion techniques. [|l9|-pl[ In this method the pair of spins 
with strongest coupling in the random chain is eliminated 
by considering the interaction with the neighboring pairs 
as a perturbation. With this procedure a new (small) 
coupling is generated between the two spins neighbor- 
ing to the eliminated pair. The iteration continues un- 
til a single pair of spins remains. As this decimation 
procedure evolves, the original initial distribution Po(K) 
of random bonds (0 < K < J) changes into a power- 
law distribution P(K) ~ A'" -1 , where a depends only 
weakly with the "cutoff" J. In this asymptotic regime, 
corresponding experimentally to sufficiently low temper- 
atures, the susceptibility behaves as x ~ T~ a , and the 
magnetic contribution to the specific heat as C ~ T 1- " 
("random singlet phase"). Since a is positive, the sus- 
ceptibility diverges as the temperature approaches zero. 
For magnetic fields much stronger than the temperature 
but much weaker than J (all in energy units), the mag- 
netization behaves as m ~ i? 1_Q (note that the actual 
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value of a is difficult to calculate within the accuracy 
of our techniques since we need to know the initial dis- 
tribution of random exchange couplings). These predic- 
tions have been verified in organic materials and also in 
MgTiOB0 3 H where it was found that a ~ 0.8. Thus, 
if indeed the low energy behavior of Zn-doped CuGeC>3 
can be mapped into a random Heisenberg model, just 
one parameter (a) is needed to fit the power-law behav- 
ior of the susceptibility, specific heat and magnetization 
as described above. This is a prediction of our theory for 
Zn-doped spin-Peierls systems in the regime of small Zn 
concentration, and it breaks down at higher concentra- 
tions when the experimentally observed transition to a 
long-range 3D Neel order takes place. 



VII. CONCLUSIONS, EXTENSION TO 
LADDERS, AND ENHANCEMENT OF SPIN 
CORRELATIONS 

Summarizing, here we have developed a simple pic- 
ture of the effect of Zn-doping on compounds with a spin 
gapped ground state. The idea is that each non-magnetic 
impurity produces an effective S=l/2 "free" spin which 
may be an actual S=l/2 electron or hole located at a 
nearest- neighbor site (in the limit of large dimerization), 
or it may correspond to many-body S=l/2 states of sev- 
ered odd-site chains (in the limit of weak dimerization). 
Similar concepts have been used in the context of other 
compounds that provide physical realizations of random 
ferro-antiferromagnetic exchange S=l/2 chains. J23] The 
interaction of the free spins is mediated by the rest of 
the spins, and thus it can be very small at low Zn- 
concentration. Such a weak coupling induces low energy 
spectral weight inside the spin gap in calculations of the 
dynamical spin structure factor. This effect occurs at all 
Zn densities, with a weight proportional to such density. 
Clear remnants of the dynamical features corresponding 
to an undoped chain are observed. On the other hand, 
the magnetic susceptibility, which reacts to low energy 
excitations, indicates the rapid suppression of the spin 
gap at very low impurity concentrations. Thus, an anal- 
ysis of INS and \ data may naively produce different 
conclusions unless we recall that INS reacts to fluctua- 
tions at all distances, while x is affected mostly by the 
long-range order in the system. 

The calculations presented in this paper have been car- 
ried out for the special case of Zn-doped CuGe03, but 
the results can be easily generalized to other system with 
spin-gaps. Of particular interest are the recent results re- 
ported in ladders, ]To[ | where a rapid collapse of the spin- 
gap was observed. To understand this result consider the 
limit where the rung exchange J± is dominant (which is 
the analog of the large <5 limit studied in Sec. II), and thus 
the ground state is made out of rung spin-singlets. [ pT| In 
this case each Cu replaced by Zn actually induces a loose 
S=l/2 spin, as recently observed numerically by Sandvik 
et al. §| (see Fig. 7). 
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FIG. 7. Schematic representation of a Zn-doped ladder. 
The open circles are spins and the full circles are impurities. 
The thick bonds are the strongest in the large J± limit. 

Their interaction should lead to small spectral weight 
in the dynamical structure factor, similar to the results 
for the SP ground state. Then, INS studies of Zn-doped 
ladders at a small Zn concentration would indicate the 
survival of weight at energies of the order of the original 
spin-gap, in regions where susceptibility measurements 
show that the gap has collapsed (which occurs rapidly 
upon doping according to the intuitive picture discussed 
in previous sections and before in Ref. ||). Recently, 
related work in the context of Zn-doped ladders have 
appeared with conclusions similar to ours. ]l3| , p5| , p6"| 

An important issue that we have not addressed in de- 
tail in this work is the behavior of the spin-spin correla- 
tions after Zn doping of CuGe03. If these correlations 
are enhanced upon such doping, then a weak interchain 
coupling may lead to the stabilization of a 3D Neel order 
as observed experimentally. 
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FIG. 8. Spin-spin correlation as a function of distance for 
a ring of 16 sites with and without impurities located at sites 
1 and 15. The parameters are J3/J1 = 0.24 and 8 = 0.03. 
The correlations are measured with respect to site 2, i.e. the 
site immediately next to an impurity. 

Actually such an unusual effect (correlation enhance- 
ment by Zn doping) has been observed in a numerical 
study of ladders fl25fl and we have also detected such an 
enhancement in one dimensional spin 1/2 chains as those 
studied here. For example, in Fig. 8 we show the spin-spin 
correlation measured on a ring of 16 sites with impurities 
located at sites 1 and 15, compared with the correlations 
obtained on a ring with the same number of sites but 
without the impurities. A remarkable enhancement of 



G 



the correlations is observed at short distances. Currently 
we are investigating the origin of this enhancement which 
seems to be generic of a variety of Zn-doped materials. 
Other systems with spin-gaps should behave similarly. 
S=l chains have been recently studied pij and weight in- 
side the Haldane gap was reported. In two dimensions, a 
similar behavior should be observed in Zn-doped cuprates 
in the temperature and density regime where a spin-gap 
is observed for the undoped case. Jl2[ 
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